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such  as  in  situations  of  severe  noise  that  are  biased  in  magnitude  and  polarity. 

The  Air  Force  Geophysics  Laboratory  was  faced  with  this  problem  in  the  analysis  of 
snow  weight /rate  data  because  of  wind -induced  extraneous  readings.  Various  automated 
analytical  schemes  were  tried  in  an  attempt  to  reproduce  the  base-line  data  but  all  failed 
to  give  approximations  within  acceptable  error  boundaries.  It  was  then  noted  that  the  base 
values  could  be  very  closely  replicated  by  a  hand -drawn,  "best  guess"  line  on  the  plotted, 
noisy-raw  data.  This  revelation  prompted  the  development  of  a  computer  process  that,  by 
mimicking  the  logic  of  human  reasoning,  can  eliminate  extraneous  readings  and  (cont. )  - 
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reconstruct  the  ba  3-line  data  to  a  very  close  approximation.  As  such,  this  computer 
decision-making  procedure  may  be  classified  as  a  form  of  artificial  intelligence  that  may 
be  applicable  to  other  analytical  routines. 

This  report  discusses  the  problems  associated  with  extraneous  or  noisy  data  and  de¬ 
scribes  the  technique  that  was  developed  to  eliminate  superfluous  readings  from  snow 
weight /rate  measurements. 
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lc.  Extraneous  Readings  Having  Symmetry  in  Magnitude 

and  Random  Polarity  Superimposed  Upon  Constant 

Base-Line  Data  2 

ld.  Extraneous  Readings  Having  Both  Random  Magnitude  and 

Polarity  Superimposed  Upon  Constant  Base-Line  Data  3 

2.  Photograph  of  Snow  Rate  Meter  8 
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3.  One  Hour  of  Snow  Weight  Data  Taken  on  10  Dec  1982 

at  Camp  Grayling,  Mich. ,  During  the  SNOW-ONE-B 
Field  Experiment 

4.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3 

5.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3 

With  the  Hand-Drawn  "Best  Guess"  Estimates  of 
the  Actual  Snow  Weights 

6.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3 

With  a  1-Min  Running-Mean  Average 

7.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3 

With  a  5-Min  Running-Mean  Average 

8.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3 

With  the  Computer-Generated  Estimates  of  the 
Actual  Snow  Weights 

9.  Flow  Diagram  Describing  the  Logic  of  the  Computer 

Decision-Making  Process 

10a.  A  Base-Line  Curve  Superimposed  on  Assumed  Computer 
Generated,  Noisy  Data  That  Varied  Randomly  in  Both 
Magnitude  and  Polarity 

10b.  The  Assumed  Noisy  Data  of  Figure  10a  With  the  Compute 
Approximation  Using  the  Noise  Elimination  Procedure 
and  a  20-Point  Running-Mean  Line 

11.  Basic  Program  of  the  Noise  Elimination  Procedure  Used 
to  Plot  the  Computer  Approximation  in  Figure  8 


A  Computer  Decision-Making  Process  for  the 
Elimination  of  Noise  from  Data 


1.  INTRODUCTION 

Data  analysis  is  often  complicated  by  the  problem  of  determining  true  values 
of  measurements  from  those  contaminated  by  extraneous  signals  normally  referred 
to  as  "noise.  "  The  most  common  method  of  alleviating  this  problem  is  to  average 
the  raw  data  over  some  period  of  time  so  that  the  resulting  mean  values  approxi¬ 
mate  the  true  readings.  This  approach  is  effective  in  many  cases  especially  when 
the  noisy  signals  are  symmetrically  superimposed  upon  a  data  base  in  both  magni¬ 
tude  and  polarity.  Data  averaging  may  also  be  effectively  employed  on  random 
noise  if  the  data  are  averaged  over  long  time  periods.  Extraneous  random  read¬ 
ings  do  not  exhibit  a  preference  in  their  deviations  from  a  data  base  and  will,  over 
a  suitably  long  time  period,  tend  to  cancel  one  antrther  to  give  accurate  mean  val¬ 
ues.  However,  if  the  extraneous  signals  are  biased  in  magnitude  and/or  polarity, 
averaging  may,  in  turn,  produce  biased  results  no  matter  how  many  points  are  in¬ 
cluded  in  the  averaging  interval.  Short-term  averaging,  on  the  other  hand,  can  be 
used  confidently  only  in  cases  where  the  superfluous  readings  are  symmetrical  in 
nature. 

The  four  plots  in  Figure  1  illustrate  the  problems  associated  with  averaging 
data  that  have  extraneous  signals  superimposed  on  base-line  data.  In  each  case. 


(Received  for  publication  26  March  1987) 


Figure  la.  Extraneous  Readings  Having  Symmetry  in 
Magnitude  and  Polarity  Superimposed  Upon  Constant 
Base-Line  Data.  A  20  point  running-mean  average 
is  indicated  by  the  straight -solid  line  and  the  base 
data  by  the  dashed  line.  The  averaged  line  corre¬ 
sponds  exactly  to  the  base-line  data 


Figure  lb.  Extraneous  Readings  Having  Random 
Magnitude  and  Symmetrical  Polarity  Superimposed 
Upon  Constant  Base-Line  Data.  A  20  point  running- 
mean  average  is  indicated  by  the  smoothed  line  and 
the  base  data  by  the  dashed  line 


i  i I'tirr  Id.  I'xlranrnus  Headings  Having  Hoth  Handom 
Magnitude  and  Polarity  Super  ini  posed  Upon  Constant 
Muse -I  me  Data.  A  20  point  running-mean  average  is 
indicated  l>v  the  smoothed  line  and  the  base  data  by  the 
dashed  lino 
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the  base  data  consists  of  200  constant  readings  as  indicated  by  the  horizontal 
dashed  lines.  Also,  in  each  case,  the  mean  value  of  the  total  200  points  (long-term 
averaging)  agrees  with  the  base  readings.  The  smoothed  lines  are  the  results  of  a 
20  point,  running-mean  average  (short-term).  Plot  la  shows  the  superfluous  read¬ 
ings  being  symmetrically  displaced  in  both  magnitude  and  polarity  from  the  base 
line.  The  short-term,  20  point  average  corresponds  exactly  to  the  base  readings. 
Plot  lb  has  symmetry  in  the  polarity  of  the  noise  (alternately  +  and  -  the  base  data) 
with  the  magnitude  of  the  deviation  selected  randomly.  The  20  point  averaged  line 
shows  small  departures  from  the  base-line  readings.  Plot  lc  shows  the  extrane¬ 
ous  signals  having  a  constant  deviation  from  the  base  data  with  a  random  polarity, 
whereas  Id  has  both  random  magnitude  and  polarity.  The  short-term  averaged 
lines  in  both  these  cases  show  a  larger  deviation  from  the  base  line  than  that  of  lb 
and  seem  to  show  that  the  effects  of  signal  polarity  can  be  more  important  than  sig¬ 
nal  magnitude. 

In  1980,  the  Cloud  Physics  Branch  of  the  Air  Force  Geophysics  Laboratory 

1  2 

designed  and  constructed  an  instrument  to  determine  the  rate  of  falling  snow.  ' 

The  operation  of  this  prototype  device  is  based  on  accurate  and  sensitive  weight 
measurements  from  an  electronic  balance  that  are  converted  into  mass  and  asso- 

3 

ciated  with  time  in  the  subsequent  calculation  of  snow  rate.  The  acquisition  of 
naturally-falling  snowflakes  requires  the  use  of  an  open  collection  container  (Fig¬ 
ure  2)  that,  in  conjunction  with  the  high -resolution  weighing  device,  makes  the  in¬ 
strument  responsive  to  the  effects  of  wind  action.  Because  of  the  ever-changing 
nature  in  the  periodicity  of  wind  gusts  and  the  variability  of  wind  velocities,  the 
raw-weight  data  display  periodic  fluctuations  of  varying  magnitude.  Several  auto¬ 
matic  reduction  schemes  were  tried  in  the  analysis  of  these  data,  but  none  could 
be  found  that  would  produce  results  within  acceptable  error  limits.  The  problem 
was  finally  resolved  with  the  development  of  a  technique  that  instructed  a  computer 
to  conduct  an  analysis  using  a  format  akin  to  human  reasoning.  Because  the  com¬ 
puter  makes  decisions  based  on  particular  characteristics  of  the  data  set  and  ad¬ 
justs  parameters  to  produce  optimum  results,  the  analytical  routine  may  be  con- 


1.  Gibbons,  L.C.,  Matthews,  A.J.,  Berthel,  R.O.,  and  Plank.  Y.G.  (1P83) 

Snow  Characterization  Instruments.  AFGI-TR -83-0063,  AD  A  13 1984. 

2.  Plank,  Y.G.,  Matthews,  A..L,  and  Berthel,  R.O.  (1983)  Instruments  used  for 

snow  characterization  in  support  of  SNOW-ONE-A  and  SNOW-ONE-B,  Pro¬ 
ceedings  of  SPIE  Technical  Symposium  East  83,  Session  'Optical  Engineer- 
ing  for  (  old  fenvironments.  "  Sub-session  Optical  Hardware  in  the  C  old,  " 
Vol.  414,  pp.  29-34. 

3.  Berthel,  R.O.,  and  Plank,  \  .G.  (1986)  Data  reduction  techniques  used  in  the 

determination  of  precipitation  rates  from  electronic  weight  measurements. 
Proceedings  of  C  orrection  of  Precipitation  Measurements,  Zurich,  Switzer- 
land,  pp.  31 -35.  ~~ 
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Figure  2.  Photograph  of  Snow  Rate  Meter.  The  outer  hinged  metal  wind 
shield  is  shown  in  the  open  position  to  expose  the  inner  plastic  shield  and 
collection  bucket 


sidered  as  a  form  of  artificial  intelligence.  Therefore,  this  technique  may  have 
beneficial  applications  in  other  analytical  and  reduction  procedures. 

The  following  sections  define  the  logic  that  led  to  the  development  of  this  tech¬ 
nique  and  describe  the  procedure  used  in  resolving  the  problem. 


2.  PROBLEM  DEFINITION 


The  data  recognition  problem  referred  to  in  the  preceding  section  became  ap¬ 
parent  from  the  analysis  of  snow  weight /rate  data.  Thus,  it  is  only  fitting  that 
this  type  data  be  used  to  define  the  problem  and  the  subsequent  solution. 

Figure  3  shows  one  hour  of  snow  weights  as  recorded  on  10  Dec  1982  during 
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Figure  3.  One  Hour  of  Snow  Weight  Data  Taken  on  10  Dec  1982  at  Camp  Grayling, 
Mich.,  During  the  SNOW-ONE-B  Field  Experiment 

4 

the  SNOW-ONE-B  field  experiment  at  Camp  Grayling,  Mich.  This  plot  consists 
of  1200  weight  readings,  one  every  3  sec,  as  measured  by  an  electronic  balance 
having  a  0.01  g  sensitivity.  The  nature  of  these  measurements,  where  the  balance 
is  weighing  the  accumulation  of  naturally-falling  snowflakes,  dictates  that  the 
weight  readings  will  remain  constant  in  the  absence  of  snowfall  and  will  increase 
during  periods  that  snow  is  falling.  The  only  allowable  deviation  from  these  rules, 
assuming  a  stable  instrument  environment  and  no  zero  drift,  is  a  decrease  in 
weight  because  of  sublimation  or  the  evaporation  of  snow  from  the  collection  con¬ 
tainer  that  is  evidenced  by  gradually  declining  weight  readings  over  long  periods  of 

4.  Bates,  R.E.,  and  Bowan,  S.L.,  Eds.  (1983)  SNOW-ONE-B  Data  Report.  U. S. 
Army  Corps  of  Engineers,  Special  Report  83-16. 


time.  Sudden  increases  or  decreases  in  weights  that  subsequently  return  to  the 
base  line  (data  spikes)  are  illogical  in  this  type  measurement  and  are  classified  as 
being  extraneous  or  data  noise.  The  base  curve  from  the  weights  of  accumulated 
snowfall  in  Figure  3  is  readily  distinguishable  as  are  the  major  positive  and  the 
minor  negative  spikes  that  are  superimposed  on  the  curve.  These  extraneous  read¬ 
ings  are  the  effects  of  wind  action  on  the  open  collection  container.  The  differences 

in  the  magnitude  of  the  spikes  above  and  below  the  data  base  line  are  attributed  to 

5  6 

a  shielding  mechanism  that  serves  to  dampen  the  negative  readings.  '  The  prob¬ 
lem,  obviously,  is  one  of  determining  the  actual  data  from  this  combination  of 
weight  and  noise  readings. 

3.  ANALYTICAL  LOGIC 

The  individual  raw -data  readings  are  more  definable  in  the  20  min  (2200  to 
2220)  plot  in  Figure  4.  This  time  segment  was  chosen  as  it  illustrates  slightly  in¬ 
creasing  weights  with  little  or  no  noise  (initial  period),  increasing  noise  and 
weights  (central  period)  and  high  noise  with  little  or  no  increase  in  weight  (latter 
period).  The  small  fluctuations  during  the  initial  quiet  period  constitute  the  nor¬ 
mal  variations  of  these  data.  By  identifying  similar  quiet  periods  between  wind 
gusts  and  drawing  a  connecting  line  through  them,  a  "best  guess"  weight  curve  may 
be  determined  as  indicated  by  the  smoothed,  hand-drawn  line  in  Figure  5.  This  is 
a  perfectly  acceptable  method  of  eliminating  the  superfluous  readings  caused  by  the 
wind  effects,  but  it  is  labor  intensive,  time-consuming,  and,  above  all,  it  is  not 
processed  in  an  automated  reproducible  manner. 

Data  smoothing  by  computer  processing  is  generally  accomplished  by  some 
form  of  data  averaging.  It  is  apparent  from  the  discussion  in  Section  1  that  the 
averaging  of  these  highly  biased,  snow-weight  data  will  not  produce  adequate  re¬ 
sults.  Figure  6  shows  the  weight  readings  from  Figure  4  with  a  superimposed 
1-min,  running-mean  line.  The  averaged  data  points  agree  favorably  during  the 
initial  3  min  of  little  or  no  wind  action,  but  show  considerable  deviation  for  the  re¬ 
mainder  of  the  time  period.  Lengthening  the  averaging  period  will  result  in  fur¬ 
ther  smoothing  as  evidenced  by  the  5 -min,  running -mean  averaged  line  in  Figure  7, 
but  the  values  will  never  match  the  true  weights  during  periods  exhibiting  strong 

5.  Berthel,  R.O.,  Plank,  V.G. ,  and  Main,  B.A.  (1983)  Analyses  of  snow  charac¬ 

terization  data  acquired  at  SNOW-ONE-A  and  B,  Snow  Symposium  III, 

CRREL,  Hanover,  N.H. 

6.  Berthel,  R.O.,  Plank,  V.  G. ,  and  Main,  B.A.  (1983)  AFGL  snow  characteri¬ 

zation  measurements  at  SNOW -ONE -B:  Preliminary  report,  SNOW  -ONE-B 
Data  Report,  U.S.  Army  Corps  of  Engineers,  CRREL  Special  Report  83-16, 
AFGL-TR-83-0174,  AD  A 130556,  pp.  197-208. 


2.3 


TIME  (Local) 

Figure  4.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3 


noisy  signals.  Long  averaging  periods  also  tend  to  blunt  or  mask  small  changes 
in  the  base  data. 

The  most  obvious  solution  to  this  problem  is  to  have  the  computer  apply  an 
analytical  process  similar  to  the  human  reasoning  used  to  construct  the  "best 
guess"  line  in  Figure  5.  To  do  this,  the  computer  has  to  be  instructed  to  find  the 
sections  of  the  raw-data  curve  where  wind  effects  were  minimal  or  not  present.  A 
straight-line  connection  of  these  quiescent  periods  should  then  approximate  the 
hand -drawn  line. 


4.  PROBLEM  RESOLUTION 

The  successful  implementation  of  this  technique  depends  entirely  upon  the 
computer  determination  of  quiescent  periods  (QPs);  therefore,  two  critical  ele¬ 
ments  of  a  QP  must  be  defined  before  this  type  of  analysis  is  conducted.  The  first 
element  is  the  minimum  time  duration  or  minimum  number  of  data  points  that  con- 
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stitute  a  QP,  and  the  second  is  the  magnitude  of  the  natural  variation  that  could  be 
experienced  in  a  situation  without  noise.  Both  are  characteristics  of  the  raw  data. 
Thus,  different  types  of  measurements  will  require  different  definitions. 

A  close  scrutiny  of  the  raw  snow-weight  data  shows  that  the  extraneous  values 
(wind  effects)  normally  exhibit  substantial  changes  in  magnitude  and  polarity  in  ad¬ 
jacent  weight  readings.  Two  noisy  readings  with  only  small  differences  between 
them  may  occur  infrequently,  but  three  adjacent  noisy  points  with  approximately 
the  same  weight  have  not  been  found.  Three  adjacent  weight  values  showing  little 
or  no  change  indicate  that  no  extraneous  readings  are  present.  The  minimum  du¬ 
ration  of  a  OP  was  therefore  defined  as  three  points  over  a  9 -sec  time  period. 

The  electronic  balance  used  in  these  measurements  had  a  resolution  of  0.  01  g. 
During  a  3 -point  period  of  no  snow  and  little  or  no  wind,  three  adjacent  readings 
would  have  three  identical  weights  or  two  of  the  same  with  the  third  showing 
+  /-  0.01  g.  Therefore,  the  natural  variation  of  the  data  allows  weight  differences 
of  at  least  +  /  -0.  01  g  within  a  QP.  This  magnitude  difference  is  referred  to  as  a 
weight  limit. 

Establishing  these  criteria  enables  the  computer  to  begin  analyzing  the  data. 
The  first  programmed  instruction  is  to  step  through  the  raw  data  and  search  for  a 
QP.  This  is  done  by  averaging  the  weight  (W)  at  a  particular  data  point  number 
(N)  with  the  preceding  weight  at  N  -  1  and  the  following  at  N  +  1  as 


W  +  W  +  W 
_  WN-1  WN  W  N  +  l 

W  = - - -  g- 
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The  mean  weight  (W)  is  then  compared  to  each  of  the  three  constituent  weights. 
If  any  of  the  three  is  found  to  be  more  than  the  designated  weight  limit  from  the 
value  of  W,  the  point  is  rejected  and  the  computer  steps  ahead  by  one  point  in  the 
raw  data  and  repeats  the  process.  When  all  three  readings  are  within  the  weight 
limit  from  W,  a  QP  has  been  found. 

At  this  time,  the  computer  must  determine  whether  this  period  is,  in  fact,  the 
next  valid  sequential  quiet  region  or  if  there  should  be  one  or  more  QPs  defined  by 
a  larger  weight  limit  between  the  QP  just  defined  and  the  last  previously  found  QP. 
This  judgment  is  necessary  so  that  the  weight  limit  may  be  adjusted  to  account  for 
increases  (or  decreases)  in  the  base-line  data.  For  example,  if  the  base  data  are 
changing  in  increments  of  0.  01  g  between  points,  a  magnitude  limit  of  0.  01  g  will 
detect  a  QP.  However,  if  the  incremental  change  were  0.  05  g,  then  a  0.  05  g  limit 
would  be  required. 

The  judgment  procedure  will  also  serve  to  adjust  the  weight  limit  to  account 
for  changes  that  may  occur  in  the  natural  variation. 


The  weight  difference  between  a  particular  QP  and  the  preceding  one  deter¬ 
mines  if  a  limit  adjustment  is  required.  If  this  differential  weight  is  larger  than 
the  magnitude  limit,  there  is  a  strong  possibility  that  the  base  data,  at  some  pe¬ 
riod  or  periods  in  the  intervening  region,  is  increasing  in  steps  larger  than  the 
designated  limit.  If  the  weight  differential  is  within  the  magnitude  limit,  that  pos¬ 
sibility  does  not  exist. 

Therefore,  when  a  quiet  period  has  been  found  using  a  particular  weight  limit, 
the  computer  determines  the  weight  difference  between  that  period  and  the  preced¬ 
ing  one.  If  the  difference  is  within  the  designated  limit  range,  the  period  is  con¬ 
sidered  valid.  If  not,  the  limit  is  increased  by  0.  01  g  and  the  search  for  a  QP  is 
begun  anew.  This  iteration  is  repeated  until  the  weight  difference  between  quies¬ 
cent  regions  is  within  the  magnitude  limit  used  to  define  the  second  QP. 

Upon  determining  the  weight  of  a  valid  quiet  period  (Wq)  at  some  data  point 
number  (Nq),  the  computer  is  then  instructed  to  construct  a  straight-line  approxi¬ 
mation  from  the  preceding  quiet  period  Wq_^  at  point  number  Nq_j.  This  opera¬ 
tion  assigns  an  interpolated  value  to  each  discarded  reading  between  the  adjacent 
QPs  and  essentially  reconstructs  the  data  curve  minus  the  noisy  points. 

The  first  step  in  this  process  is  the  determination  of  the  incremental  weight 
difference  (AW)  that  is  to  be  applied  between  data  points  by 

(2) 
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The  individual  intervening  points  are  then  assigned  a  weight  (W^)  as 


N, 
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WN  =)WQ-1  +AW  <N  -NQ.^g. 


(3) 
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(In  this  description,  the  weight  readings  are  associated  with  data  point  numbers 
since  the  3-sec  time  increments  between  readings  may  possibly  cause  confusion  if 
the  weights  were  based  on  a  time  reference. ) 

When  the  raw  data  from  the  2200  to  2220  time  period  were  subjected  to  the 
above  analysis,  a  computer-generated  approximation  of  the  true  data  curve  was 
constructed  as  shown  by  the  superimposed  line  in  Figure  8.  The  difference  in  this 
line  and  the  "best  guess"  line  of  Figure  5  are  minor  and  few.  One  can  only  specu¬ 
late  as  to  which  approximation  is  more  exact.  The  hand -drawn  line  tends  to  be 
smoother  as  some  slight  dips  and  rises  can  be  detected  in  the  computer  version. 
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Figure  8.  A  20-Min  Segment  of  the  Noisy  Data  From  Figure  3  With  the  Computer- 
Generated  Estimates  of  the  Actual  Snow  Weights 
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The  computer  decision-making  process  discussed  herein  and  described  in  the 
flow  diagram  of  Figure  9  was  purposely  designed  to  eliminate  wind-induced  extra¬ 
neous  readings  from  weight  measurement'  used  in  the  determination  of  snow  rate.  , 
It  must  be  emphasized  that  the  application  of  this  procedure  to  other  type  measure¬ 
ments  imposes  two  requirements  on  the  raw  data.  First,  the  data  must  be  re¬ 
corded  or  stored  since  iteration  may  be  necessary,  and,  second,  the  data  must 
periodically  approximate  the  base-line  i  eadings  or,  in  other  words,  must  have 
quiet  or  minimal  noisy  periods. 

In  addition,  it  may  be  advantageous  to  designate  a  particular  point  limit  in  the 
search  for  QPs.  In  most  cases,  the  OP  validity  check,  where  the  differences  in 
QP  values  are  required  to  conform  to  the  magnitude  limit,  will  act  to  restrict  the 
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Figure  9.  Flow  Diagram  Describing  the  Logic  of 
the  Computer  Decision-Making  Process.  Defini¬ 
tions  of  terms  are  found  in  Section  4 

number  of  points  in  the  search  for  a  QP  because  of  changing  base-line  data.  How¬ 
ever,  in  certain  situations  of  heavy  noise  and  a  slowly  changing  base  line,  a  data 
set  may  conceivably  be  described  by  one  straight-line  extrapolation  over  a  large 
number  of  points.  This  circumstance  is  alleviated  by  designating  a  maximum  time 
duration  or  maximum  number  of  points  in  the  search  for  a  QP.  This  designation 
is  referred  to  as  a  "search  limit.  " 

Upon  reaching  the  number  of  points  in  the  search  limit,  it  is  assumed  that  any 
and  all  intervening  quiet  periods  have  larger  deviations  than  that  of  the  current 
magnitude  limit.  This  limit  is  then  increased  and  the  search  restarted. 

Snow  weight  data  generally  show  quiet  regions  occurring  periodically  in  'v  1  to 
^  4  min  intervals.  A  5-min  or  100-point  search  limit  was  included  in  the  snow  rate 
analysis  program  and  in  the  noise  elimination  process  used  to  produce  the  comput¬ 
er  generated  approximation  in  Figure  8.  However,  the  limit  was  not  required  in 
this  20-min  segment  as  the  maximum  search  period  encompassed  85  points  (4.25 
min). 

The  noise  elimination  process  has  not  been  applied  to  measurements  other 
than  snow  weights,  although  it  has  proved  successful  in  tests  using  computer¬ 
generated  assumed  data.  For  example.  Figure  10a  shows  200  points  of  base-line 
data  in  the  form  of  a  partial  sine  wave  with  superimposed  assumed  noisy  "raw 
data. "  The  extraneous  readings  are  random  in  both  magnitude  and  polarity.  Also, 
the  QPs  were  randomly  dispersed  at  intervals  of  <  30  points,  thus  imposing  a  30- 
point  search  limit.  The  computer  approximation  of  the  base  data  is  shown  in  Fig¬ 
ure  10b  along  with  a  20-point  running-mean  average  of  the  assumed  raw  data. 
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Figure  10a.  A  Base-Line  Curve  Superimposed  on  Assumed 
Computer-Generated,  Noisy  Data  That  Varied  Randomly  in 
Both  Magnitude  and  Polarity 


Figure  10b.  The  Assumed  Noisy  Data  of  Figure  10a  With  the 
Computer  Approximation  Using  the  Noise  Elimination  Proce¬ 
dure  and  a  20-Point  Running-Mean  Line 
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The  procedure  described  in  this  report  is  straightforward  in  logic  and  appli¬ 
cation.  Because  of  this,  it  may  prove  useful  in  the  analysis  of  other  types  of 
measurements.  The  entire  process  requires  a  relatively  simple,  short  computer 
routine  as  demonstrated  by  the  27-line  basic  program  in  Figure  11.  This  program 
was  used  to  determine  the  values  of  the  computer-generated  approximation  shown 
plotted  in  Figure  8. 


D(n)  =  Data  (raw)  array 

X(n)  =  Point  number  array  for  plotting 
V(n)  =  Weight  array  for  plotting 

A  =  Average  of  3  points 

L  =  Magnitude  limit 

N  Number  of  plotting  points 
P  =  Point  number  at  start  of  OP  search 
W  =  W  eight  at  preceding  OP 
W  1  -  Delta  weight  between  adjacent  OPs 


100 

N  =  1 

1  10 

P-1 

120 

L  - .  0 1 

(Iteration) 

130 

For  I-P+l  to  399 

140 

A  =  (D(I-1  )  +  D(l)  +  D(I-l  ))/3 

(Eq.  1) 

150 

For  3=1=1  to  I+1 

160 

If  A  -D(J )  >  L  then  350 

(check  search  limit) 

170 

If  D(J )-A  >  L  then  350 

(check  search  limit) 

180 

Next  J 

(QP  found) 

190 

If  N  >  I  then  230 

(validity  check) 

(initial  plotting  point) 

200 

X  ( N )  =  I 

210 

V(N)  =  A 

220 

Goto  320  (continue  sampling) 

(Validity  check) 

230 

If  A  -W  1 .  then  260 

(valid  OP) 

(Increasing  magnitude  limit) 

240 

LI. +  .01 

250 

Goto  130  (iteration) 

(Valid  QP) 

260 

W  1=(A-W)/(I-P>  (Eq.  2) 

(Interpolation) 

270 

For  J  =  P  +  1  to  I 

280 

N  =N  +  1 

Figure  11.  Basic  Program  of  the  Noise  Elimination 
Procedure  Used  to  Plot  the  Computer  Approximation 
in  Figure  8 
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I  igure  11  (cont.  ).  Basic  Program  of  the  Noise  Elimination 
Procedure  l  set!  to  Plot  the  Computer  Approximation  in  Fig¬ 
ure  1! 
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